Lecture 13: A low-rank factorization framework for building scalable algebraic solvers and preconditioners

Factorization based preconditioning algorithms, most notably incomplete LU (ILU) factorization, have been shown to be robust and applicable to wide ranges of problems. However, traditional ILU algorithms are not amenable to scalable implementation. In recent years, we have seen a lot of investigations using low-rank compression techniques to build approximate factorizations.A key to achieving lower complexity is the use of hierarchical matrix algebra, stemming from the H-matrix research. In addition, the multilevel algorithm paradigm provides a good vehicle for a scalable implementation. The goal of this lecture is to give an overview of the various hierarchical matrix formats, such as hierarchically semi-separable matrix (HSS), hierarchically off-diagonal low-rank matrix (HODLR) and Butterfly matrix, and explain the algorithm differences and approximation quality. We will illustrate many practical issues of these algorithms using our parallel libraries STRUMPACK and ButterflyPACK, and demonstrate their effectiveness and scalability while solving the very challenging problems, such as high frequency wave equations

The YH database: the first Asian diploid genome database

The YH database is a server that allows the user to easily browse and download data from the first Asian diploid genome. The aim of this platform is to facilitate the study of this Asian genome and to enable improved organization and presentation large-scale personal genome data. Powered by GBrowse, we illustrate here the genome sequences, SNPs, and sequencing reads in the MapView. The relationships between phenotype and genotype can be searched by location, dbSNP ID, HGMD ID, gene symbol and disease name. A BLAST web service is also provided for the purpose of aligning query sequence against YH genome consensus. The YH database is currently one of the three personal genome database, organizing the original data and analysis results in a user-friendly interface, which is an endeavor to achieve fundamental goals for establishing personal medicine. The database is available at http://yh.genomics.org.cn

Differential Developmental Deficits in Retinal Function in the Absence of either Protein Tyrosine Sulfotransferase-1 or -2

To investigate the role(s) of protein-tyrosine sulfation in the retina and to determine the differential role(s) of tyrosylprotein sulfotransferases (TPST) 1 and 2 in vision, retinal function and structure were examined in mice lacking TPST-1 or TPST-2. Despite the normal histologic retinal appearance in both Tpst1−/− and Tpst2−/− mice, retinal function was compromised during early development. However, Tpst1−/− retinas became electrophysiologically normal by postnatal day 90 while Tpst2−/− mice did not functionally normalize with age. Ultrastructurally, the absence of TPST-1 or TPST-2 caused minor reductions in neuronal plexus. These results demonstrate the functional importance of protein-tyrosine sulfation for proper development of the retina and suggest that the different phenotypes resulting from elimination of either TPST-1 or -2 may reflect differential expression patterns or levels of the enzymes. Furthermore, single knock-out mice of either TPST-1 or -2 did not phenocopy mice with double-knockout of both TPSTs, suggesting that the functions of the TPSTs are at least partially redundant, which points to the functional importance of these enzymes in the retina

Local alignment of two-base encoded DNA sequence

<p>Abstract</p> <p>Background</p> <p>DNA sequence comparison is based on optimal local alignment of two sequences using a similarity score. However, some new DNA sequencing technologies do not directly measure the base sequence, but rather an encoded form, such as the two-base encoding considered here. In order to compare such data to a reference sequence, the data must be decoded into sequence. The decoding is deterministic, but the possibility of measurement errors requires searching among all possible error modes and resulting alignments to achieve an optimal balance of fewer errors versus greater sequence similarity.</p> <p>Results</p> <p>We present an extension of the standard dynamic programming method for local alignment, which simultaneously decodes the data and performs the alignment, maximizing a similarity score based on a weighted combination of errors and edits, and allowing an affine gap penalty. We also present simulations that demonstrate the performance characteristics of our two base encoded alignment method and contrast those with standard DNA sequence alignment under the same conditions.</p> <p>Conclusion</p> <p>The new local alignment algorithm for two-base encoded data has substantial power to properly detect and correct measurement errors while identifying underlying sequence variants, and facilitating genome re-sequencing efforts based on this form of sequence data.</p

De novo mutations in SMCHD1 cause Bosma arhinia microphthalmia syndrome and abrogate nasal development

Bosma arhinia microphthalmia syndrome (BAMS) is an extremely rare and striking condition characterized by complete absence of the nose with or without ocular defects. We report here that missense mutations in the epigenetic regulator SMCHD1 mapping to the extended ATPase domain of the encoded protein cause BAMS in all 14 cases studied. All mutations were de novo where parental DNA was available. Biochemical tests and in vivo assays in Xenopus laevis embryos suggest that these mutations may behave as gain-of-function alleles. This finding is in contrast to the loss-of-function mutations in SMCHD1 that have been associated with facioscapulohumeral muscular dystrophy (FSHD) type 2. Our results establish SMCHD1 as a key player in nasal development and provide biochemical insight into its enzymatic function that may be exploited for development of therapeutics for FSHD

SeqGene: a comprehensive software solution for mining exome- and transcriptome- sequencing data

Abstract Background The popularity of massively parallel exome and transcriptome sequencing projects demands new data mining tools with a comprehensive set of features to support a wide range of analysis tasks. Results SeqGene, a new data mining tool, supports mutation detection and annotation, dbSNP and 1000 Genome data integration, RNA-Seq expression quantification, mutation and coverage visualization, allele specific expression (ASE), differentially expressed genes (DEGs) identification, copy number variation (CNV) analysis, and gene expression quantitative trait loci (eQTLs) detection. We also developed novel methods for testing the association between SNP and expression and identifying genotype-controlled DEGs. We showed that the results generated from SeqGene compares favourably to other existing methods in our case studies. Conclusion SeqGene is designed as a general-purpose software package. It supports both paired-end reads and single reads generated on most sequencing platforms; it runs on all major types of computers; it supports arbitrary genome assemblies for arbitrary organisms; and it scales well to support both large and small scale sequencing projects. The software homepage is http://seqgene.sourceforge.net.</p

A Bi-Functional Anti-Thrombosis Protein Containing Both Direct-Acting Fibrin(ogen)olytic and Plasminogen-Activating Activities

Direct-acting fibrin(ogen)olytic agents such as plasmin have been proved to contain effective and safety thrombolytic potential. Unfortunately, plasmin is ineffective when administered by the intravenous route because it was neutralized by plasma antiplasmin. Direct-acting fibrin(ogen)olytic agents with resistance against antiplasmin will brighten the prospect of anti-thrombosis. As reported in ‘Compendium of Materia Medica’, the insect of Eupolyphaga sinensis Walker has been used as traditional anti-thrombosis medicine without bleeding risk for several hundreds years. Currently, we have identified a fibrin(ogen)olytic protein (Eupolytin1) containing both fibrin(ogen)olytic and plasminogen-activating (PA) activities from the beetle, E. sinensis. Objectives: To investigate the role of native and recombinant eupolytin1 in fibrin(ogen)olytic and plasminogen-activating processes. Methods and Results: Using thrombus animal model, eupolytin1 was proved to contain strong and rapid thrombolytic ability and safety in vivo, which are better than that of urokinase. Most importantly, no bleeding complications were appeared even the intravenous dose up to 0.12 µmol/kg body weight (3 times of tested dose which could completely lyse experimental thrombi) in rabbits. It is the first report of thrombolytic agents containing both direct-acting fibrin(ogen)olytic and plasminogen-activating activities. Conclusions: The study identified novel thrombolytic agent with prospecting clinical potential because of its bi-functional merits containing both plasmin- and PA-like activities and unique pharmacological kinetics in vivo

Specialized dynamical properties of promiscuous residues revealed by simulated conformational ensembles

The ability to interact with different partners is one of the most important features in proteins. Proteins that bind a large number of partners (hubs) have been often associated with intrinsic disorder. However, many examples exist of hubs with an ordered structure, and evidence of a general mechanism promoting promiscuity in ordered proteins is still elusive. An intriguing hypothesis is that promiscuous binding sites have specific dynamical properties, distinct from the rest of the interface and pre-existing in the protein isolated state. Here, we present the first comprehensive study of the intrinsic dynamics of promiscuous residues in a large protein data set. Different computational methods, from coarse-grained elastic models to geometry-based sampling methods and to full-atom Molecular Dynamics simulations, were used to generate conformational ensembles for the isolated proteins. The flexibility and dynamic correlations of interface residues with a different degree of binding promiscuity were calculated and compared considering side chain and backbone motions, the latter both on a local and on a global scale. The study revealed that (a) promiscuous residues tend to be more flexible than nonpromiscuous ones, (b) this additional flexibility has a higher degree of organization, and (c) evolutionary conservation and binding promiscuity have opposite effects on intrinsic dynamics. Findings on simulated ensembles were also validated on ensembles of experimental structures extracted from the Protein Data Bank (PDB). Additionally, the low occurrence of single nucleotide polymorphisms observed for promiscuous residues indicated a tendency to preserve binding diversity at these positions. A case study on two ubiquitin-like proteins exemplifies how binding promiscuity in evolutionary related proteins can be modulated by the fine-tuning of the interface dynamics. The interplay between promiscuity and flexibility highlighted here can inspire new directions in protein-protein interaction prediction and design methods. © 2013 American Chemical Society